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ABSTRACT 

Aims. We investigate the large scale metallicity distribution in the Galactic bulge, using a large spatial coverage, in order to constrain 
the bulge formation scenario. 

Methods. We use the VISTA variables in the Via Lactea (VVV) survey data and 2MASS photometry, covering 320 sqdeg of the 
Galactic bulge, to derive photometric metallicities by interpolating of the (J - Ks) colors of individual Red Giant Branch stars based 
on a set of globular cluster ridge lines. We then use this information to construct the first global metallicity map of the bulge with a 
resolution of 30' x 45'. 

Results. The metallicity map of the bulge revealed a clear vertical metallicity gradient of ~ 0.04dex/deg (~ 0.28dex/kpc), with metal- 
rich stars ([Fe/H]~0) dominating the inner bulge in regions closer to the galactic plane (\b\ < 5). At larger scale heights, the mean 
metallicity of the bulge population becomes significantly more metal-poor. 

Conclusions. This fits in the scenario of a boxy-bulge originated from the vertical inestability of the Galactic bar, formed early via 
secular evolution of a two component stellar disk. Older, metal-poor stars dominate at higher scale heights due to the non-mixed orbits 
from the originally hotter thick disk stars. 
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• , 1. Introduction collisions and mergers, stars in classical bulges are expected to 

£^ . have formed mostly from the starburst in the early times, with 

O • In the hierarchical scenario of galaxy formation, the early mass little or no star formation since then. For this reason, classi- 

CO ; assembly of galaxies will be dominated by merger processes cal bulges are expec ted to be composed predominantly of old 

. which favors the formation of spheroids (|Toomre||1977|). This stars ( ^ 10 G vrl (e.g. lRenzinilll998l they have enhanced abun- 

■ violent and relatively fast process produces the present day el- dance ratios dPeletier et alJ ll999r. iMoorth v & Holtzmanl I2006I 

k> : liptical S alaxies - However, when the spheroids are accompa- |MacArthur et al.l2008l) and often show metallicity gradients, go- 

k£ . nied by the (re-)growth of stellar disk flKauffmann et al. || 1999t ing from me tal-poor stars in the outskirts to a more metal-rich 

g ; ISpnngel & Hernqmst||2005|), they receive the name of classical population towards the center, (see for example Uablonka et all 

■ " " bul 8 es - l2007h . 

Classical bulges are expected to share the same observational Given that the classical bulge formation scenario is con- 
properties of elliptical galaxies. Their light profil es are well fit- nec ted to the rate of merger events, this is expected to be 
ted by a Vaucouleurs-like profile with n~4 (e.g. |Gadotti||2009l) the do minant process in the evolut ion of galaxies at early 
and their sta rs show hot d ynamics dominated by velocity dis- times dKormendv & Kennicuttl l2004h . Later on, secular evo- 
persion (e.g.|Gadotti|[201 1[). An important point is that the rota- i ution proce sses would become important, re-arranging stel- 
tion of classical bulges is observed to be non-cylindrical, there- ^ material from a settled galactic disk. Disks instabil- 
fore showing a decrease i n the rotation velocity at increasing itieSi such as the presence of spiral arms, are expected 
height from the plane (e.g. |Emsellem et aD|2004|). Since classi- to induce the f ormation of bars in the inner regions of 
cal bulges are formed violently, in fast and early processes of a galactic disk dCombes & Sanders! 1 198 U IWeinberg] fl985; 

iDebattista & Sellwoodlll998r. lAthanassoula & Misiriotisll2002l) . 

Send offprint requests to: O. A. Gonzalez Indeed, galactic bars are a very common phenom ena, found in ~ 

* Based on observations taken within the ESO VISTA Public Survey 60% of spiral galaxies in the local Universe (e.g. lEskridge etail 

VVV, Programme ID 179.B-2002 l2000h . 
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Although disks are thin galactic components, bars might be- 
come thicker because of the vertical resonances and buckle off 
the plane of the galaxy. For this reason, galaxies, when ob- 
served edge-on, often show a central component which swells 
out of the disk with a boxy, peanut, or even an X-shape 
(e.g. iBureau et ail l2006t iPatsis et all 120021) . These structures 
have been seen in numerous simulations, as the result of the 
developmen t of vertical instabi l ities from the, initially thin, 
stellar bar (lAfh anassoul al l2005t IDebattista & Williams! 2004; 
IDebattista et al.ll2006tlMarfinez-Valpuesta et al.ll2006l) . 

The dynamical and structural properties of boxy-bulges are 
well constrained from models and observations. However, the 
properties of their stellar populations are less clear than in the 
case of classical bulges. In fact, a variety of stellar popula- 
tions are observe d within bars (e.g. iGadotti & de Souzal 12006; 
Perez et al. 2009). In practice, since boxy-bulges are bars thick- 
ened by vertical instabilities, the null hypothesis would be that 
they should hold similar ages as their parent disks. Similar case 
would be expecte d for the metalli city and alpha elements abun- 
dances. However. iFreemanl (12008) suggested that age and abun- 
dance gradients could be present in boxy bulges. Older, more 
metal-poor stars would have more time to be scattered to larger 
distances from the plane by the buckling process, therefore es- 
tablishing a negative metallicity and positive age gradient. The 
presence and strength of such gradients would then depend on 
the chemical enrichment history of the parent disk material, and 
the timescale of the buckling instability to form the boxy-bulge. 

The brightness profiles of boxy-bulges are less centrally 
concentrated than those of classical bulges, and are bet- 
ter fitted by light profiles with a lower s ersic index (e.g. 
IKormendv & Bruzual A.I \l9T& iGadottil l2009h . What is proba- 
bly the most characteristic property of boxy-bulges its that they 
present cylindrical rotation patterns. This means that the rota- 
tion velocity within a boxy-bulg e does not change with height 
above the galactic plane (e.g. IKormendv & Illingworth 1982; 
Com bes et al]1990tfAlhanassoula' & Misiriotis 2002; Shen et al.1 



2010). 



The presence of a central bar has important consequences for 
the subsequent evolution of the galaxy, redistributing the angular 
momentum and matter in the disk. In particular, it causes the loss 
of angular momentum of the molecular gas settled in the rotat- 
ing disk, inducing its infall process into the inner regions of the 
galaxy. The gas, slowly fueling into the center, produces a slow 
star formation process that results in the formation of a disk-like 
component in these inner regions. We refer to these structures 
as pseudo-bulges. The brightness of pseudo-bulges can be fit- 
ted with a very disk-like, exponential light profile and they are 
observed as flat, rotation supported structures. Also, as a result 
of a very slow process of gas fueling into the center because of 
the presence of the bar, the star formation is expected to be slow 
and continuous. A direct consequence is that pseudo-bulges host 
young stellar populations and are dominated by metal-rich stars 
with solar alpha element abundances. 

As described before, boxy-bulges are thickened bars, with 
bars being in fact disk phenomena resulting from internal pro- 
cesses of secular evolution. The absence of an external factor in 
the formation process of boxy-bulges and the fact that they are 
both formed our of originally disk material, as opposite to the 
classical bulge scenario, led different authors to refer to them 
as pseudo-bulges. This was also qualit atively supported by the 
pseud o-bulge definition in the review of IKormendv & Kennicutll 
(2004), where they describe: '..if the component in question is 
very E-like we call it bulge, if its more disk-like, we call it 
pseudo-bulge' . The morphology of thickened bars is certainly 



not E-like (Burea u et al . 2006). However, to avoid confusion 
during our analysis, we refer to them as boxy-bulges and not 
pseudo-bulges. 

The above scenarios for bulge formation and their properties, 
mostly based on the study of bulges in nearby galaxies, has been 
recently challenged by the observation of galaxies at redshift 
~ 2, i.e., at a lookback time comparable to the ages seen in clas- 
sical bulges. At this redshift, disk galaxies are observed to hold 
large, rotating disks with a very high gas fraction and velocity 
dispe r sion, significatly different from local spirals (Genze l et al.l 
2006, l2008| (Forster Schreiber et all 120091 iTacconi et al.ll2010t 
Dad di et al J2010I) . Massive gas clumps are forming in these gas- 
rich galaxies via disk intestabilities, and then interact and likely 
merge in the center of the disk. This violent process naturally 
leads to the rapid formation of a bulge over timescales of a few 
10 8 yr, much shorter than that typically ascribed to secular insta- 
bilities in local disk galaxies dlmmeh et alJl2004l : ICarollo et al.l 
120071; lElmegreen et al.l 12008). and therefore results in an old, 
ff-element enhanced, p resent day stellar populations dRenzinil 
120091: iPeng etalj 12010) of both bulge and disk formed at the 
early cosmic epoch. 

With the exception of the Milky Way and M31 
dJablonka & Sarajedinil 120051) . the observational studies of 
the central components of galaxies, based on which most 
of the above formation scenarios are derived, are strongly 
challenged by the impossibility to resolve individual stars due to 
crowding and blending in high stellar density regions of distant 
galaxies. In this context, the Milky Way bulge, one of the major 
components of the Galaxy, can be studied at a unique level 
of detail thanks to its resolved stellar populations which hold 
the imprints of how our Galaxy formed and evolved. Thus, a 
better understanding of the Milky Way bulge will provide solid 
constraints to interpret observations of extragalactic bulges as 
well as for detailed galaxy formation models. 

The place of the Milky Way within these different bulge for- 
mation scenarios is far from being understood. The Galactic 
bulge st ructure is observed to be dominated by a stel- 
lar bar (IStanek et al.l 1 19971) oriented ~ 25 d eg with re- 
spect to the Sun-Galactic center line of sight (Stanek et al. 
1 19971: lL6pez-Corredoira et al.l 120051: [Rattenburv et al.l 120071 
bonzalez et al.l l2012t iNataf et al.l 120121) . 
the bulge stars is 10 Gyr or older, based on t he deep pho- 



tometric studies in few low exti nction windows (IZoccali et al.l 
120031; IClarkson etaT1l201ll) . but iBensbv et all (120121) found a 
non-negligible fraction of younger stars based on a microlensing 
study. A vertical met allicity gradient has been observed along 
the bulge minor axis dZoccali et al.l 120081: iJohnson et al.l 1201 ll) 



that most likely flattens in the inner regions dRich et all 2007b 



In addition to the population of stars associated with a metal- 
rich, alpha-poor stellar bar, there are hints for a distinct compo- 
nent dominated by more metal-poor, alpha enriched stars show- 
ing hot kinematics wh i ch resemble those of a classica l bulge 
(Bab usiaux et"aT)l20Tot Iffill et alJ[2QTTt iGonzalez et all 12011 at 
lUttenthaler et al.l 120121) . Furthermore, significant evidence is 
now available suggest ing that the Galactic bulge i s X-shaped 
in the outer regions dMc William & Zoccalil 1201 Ot ISaito et ail 
1201 ll) . as expected from buckling instability processes of the bar, 
although the metal- poor stars do no t seem to follow the same 
spacial distribution (Nes s et al.ll2.Q12b . 

Such complex set of properties for the Galactic bulge are 
the result of the detailed study of discrete fields, mostly con- 
centrated along the minor axis, especially for the case of the 
chemical abundances. A larger, more general mapping of these 
properties is the necessary bridge between the detailed studies of 
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the Galactic bulge and those of external bulges, where integrated 
light studies can only provide us with the global information. 

In iGonzalez et alJ (1201 lbl) (Paper I) we have shown a pro- 
cedure to build multiband near-IR catalogs using the VISTA 
Variables in the Via Lactea (VVV) public survey data. We have 
then used these to derive a complete, high-resolution extinction 
map of the Galactic bulge, based on the color information of 
the red c lump (RC) giant stars . The extinction map was pre- 
sented in Gon zalez et alJ (120121) (Paper II) together with a dis- 
cussion of the important implications of its use in bulge studies. 
One of the improvements is the availability of a de-reddened RC 
feature in the luminosity function which allows us to trace the 
mean distance of stars towards different regions of the bulge, the 
so-called red clump method. The mean distance to each line of 
sight allows us then to place our photometric catalogs on the ab- 
solute plane and to derive photometric metallicity distributions. 
In Paper I we have tested the method by comparing these photo- 
metric metallicity distributions for regions along the bulge minor 
axis, where accurate spectroscopic metallicities are available, 
finding a remarkable agreement between both methods. In this 
article, the last of the series, we extend the derivation of photo- 
metric metallicities to the rest of the region covered by the VVV 
survey. We provide for the first time the global view of the metal 
content of bulge stars and its gradients as if we were looking at 
the Milky Way as an external galaxy. 



3. Extinction and mean distances 

Applying the extinction correction derived for the whole bulge 
on a spatial resolution of 2' (\b\ < 4), 4' (4 < \ b\ < 6), or 6' 
(b > 6 ') (Paper II) and using the extinction law of ICardelli et al.l 
d 1989!) . the de-reddened luminosity function is derived for each 
line of sight in the bulge. The distances to these lines of sight are 
then derived by applying the equation: 

(m - M) = K S (RC) - M Ks (RC) - A Ks + AM Ks ( 1 ) 

Mk, = -1.55 is adopted for the R C luminosity for 10 Gyr 
old population with a solar metallicity (Salaris & Girardi 2002; 
Pietrinfe mi et al.ll2004t Ivan H elshoech t & Groenewegenl l2007). 
which is an avera ge value appropriate for the bulk of the Milky 
Way bulge stars dZoccali et al.l 1200 31). However, the intrinsic 
magnitude of the RC (Mk s ) is known to have a certain depen- 
dence on metallicity. Although this dependence is neglected by 
most studies, a large metallicity gradient present in the bulge 
could bias our own metallicity calculation if a single Mk s value 
is adopted for the whole bulge. In order to take into the account 
the (expected) presence of the metallicity variations we apply 
the metallicity correction term (AM^ s ). As a consequence, an 
iterative process, described in the next section, is adopted to de- 
rive the fully consistent distance and metallicity for each line of 
sight. 



2. The data 



We use the VVV public DR1 catalogues lSaito et al.l (l2012h . For 
a complete description of the survey and the observation strategy 
we refer the reader to iMinniti et al.l (120101) . A detailed descrip- 
tion of the combination of the tile-based source lists and their 
calibration is given in Paper I. Methods used in Paper I and Paper 

11 are based on the analysis of the RC properties, which are re- 
stricted to a magnitude range well covered by the VVV survey 
photometry, across ~ 320 sqdeg of the Galactic bulge, spanning 
from -10 < / < 10 and -10 < b < +5. For this final part of the 
project, we require upper RGB photometry to derive photomet- 
ric metallicities by color interpolation be tween globular clusters 
ridge lines. As d escribed in other articles (Gonza lez et al.l2 01 lb; 
Saito et al1l2012l) . VVV photometry suffers from saturation for 
stars brighter than K s — 12 which corresponds to the upper RGB 
region in the bulge CMDs required for the color interpolation. 
For this reason, in Paper I, we have constructed VVV multiband 
catalogs and calibrated them into the 2MASS system, in order 
to later replace the saturated range of magnitudes with those of 
2MASS photometry. 

The VVV-2MASS calibration procedure is based on the 
comparison of stars flagged in 2MASS with high photomet- 
ric quality in all three bands, in the magnitude range between 

12 < K s < 13. This magnitude range avoids the saturated stars 
in VVV while being bright enough to ensure a sufficient num- 
ber of 2MASS sources with accurate photometry for the calibra- 
tion. The zero-point was calculated and applied individually to 
each VVV tile catalog. Finally stars with K s < 12 mag were re- 
moved from the VVV catalogue and were replaced with 2MASS 
sources. 

Since 2MASS is severly affected by crowding in the regions 
located close to the galactic plane, we restrict our analysis for 
latitudes \b\ > 2. This restriction aids also to avoid the high ex- 
tinction regions where some differential reddening, in variation 
scales smaller than 2', the resolution of the bulge extinction map 
from Paper II. 



4. Photometric metallicities 

The method to derive photometric metallicities in this work is 
based on the interpolation of the dereddened [J - K s )q colors 
of upper RGB stars between ridge lines of globular clusters with 
well-know metallicities derived from spectroscopy. In Paper I we 
have compared the photometric metallicity distribution along the 
bulge minor axis, to that of high resolution spectroscopy for the 
same regions. We observed a very good agreement between both 
results, and in particular, the photometric results reproduced the 
minor axis metallicity gradient that has been confirmed in nu- 
merous specotroscopic studies. 

The complete method can be summarized as follows: 

1 . The K So luminosity functions are built for each VVV tile and 
the mean magnitude of the RC is calculated by fitting the 
equation: 



N(K So ) 



a + bK Sa + cK; o + 



RC 



VRC 



V2^ 



exp 



2(r lc 



(2) 



2. Each tile is divided in four subfields of 30'x45'. The distance 
modulus is then derived for each subfield by (K s o - Mk s )rc 
using equation 1 . In the first step, AMj;, = 0, but this value 
changes in subsequent iterations. 

3. The (J - K s ) vs. M Ks CMD for each 30' x 45' sub-field is 
used to select stars with M Ks > -4.5 and (J - K s ) > 0.65, 
selecting the RGB region above the AGB bump that is more 
sensitive to metallicity, and avoiding most of the foreground 
disk contamination. An upper magnitude cut is also applied 
at the Mk, corresponding t o the RGB tip meas urements of 
bulge globular clusters from lValenti et al.l (12007). 

4. The color of each star is then interpolated among the em- 
pirical templates of Galactic halo and bulge globular clus- 
ters (M92, M55 , NGC 6752, NGC 362, M 69, NGC 6440, 
NGC 6528 and NGC 679 1 ), with well known metallicity, se- 
lected from the sample of lValenti et al.l (12007). An individual 
metallicity value is obtained for the stars in each subfield by 
this color interpolation. 
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is the result of a central concentration of metal-rich stars in the 
region -5° < b < 5° and -7° < I < 7°. This concentration of 
near-solar metallicity stars seems to follow a radial gradient of 
~ 0.04 dex/deg, going from more metal-poor stars in the outer 
bulge to more metal-rich stars towards the inner regions. 




^^HM -11.55 ll 



-2.0 -1.5 



-0.5 0.0 
[Fe/H] 



Galactic Longitude (1) 

Fig. 2. Map of the mean values of the metallicity distributions for 
the Galactic bulge covered by the V V V survey using the Cardelli 
extinction law. 



As described in Paper I, the metallicity map will be publicly 
available to the community via the BEAM (Bulge Extinction 
And Metallicity) calculator webpag^J 



Fig. 1. Color-magnitude diagram of the Galactic bulge in 
Baade's window (30' x 30'), in the selected upper RGB region, 
used for the interpolation routine. Stars are color coded acording 
to their individual metallicities. Lower panel shows the MDF for 
the same region, compared to that from high-resolution spec- 
troscopy from Hill et al. (201 1) shown in red. 



5. The mean metallicity of the subfield is then obtained and the 
AMk s va lue is calculated, using the new mean metallicity, 
based on Salaris & Girardi (2002) equation: 

AM Ks = Q.275[Fe/H] + 1.496 - M Ks (RC), (3) 

and applied to equation 1 . 

Steps 2-5 are iterated to obtain a final mean metallicity, that 
corresponds to a single Mk,. An example of the CMD in the 
Baade's window, color coded by the metallicity of the individual 
stars, is shown in the upper panel of Figure 1 . The lower panel 
compares the photometric metallicities (hashed histog ram) with 
the hig h resolution spectroscopic measurements from iHill et al.l 
(201 1) presented in red. 

Having proven the reliability of the mean metallicity mea- 
surements from the VVV+2MASS photometry (Paper I) we can 
now present in Fig. [2] the first global map of mean metallicity 
in bins of 30' x 45'. This map is restricted to \b\ > 3° as the 
upper RGB in our dataset corresponds to 2MASS photometry 
which is limited by crowding in the inner regions. This is the 
first global view of the mean metallicity variations across the 
Milky Way bulge. The gradient along the minor axis, observed 
from high-resolution studies, is clearly seen in our map and it 



5. Sources of error 

Clearly, the most accurate abundance measurements are 
obtained using high-resolution spectroscopic methods. 
Unfortunately, such studies are very expensive in terms of 
telescope time, allowing only a few regions to be studied and 
thus not allowing us to obtain the large scale picture of the bulge 
chemical abundance properties. The accuracy of the photometric 
abundance derivations depend strongly on the knowledge of 
other properties affecting the stellar population under study such 
as the density distribution in the corresponding line of sight 
and the interstellar extinction. Although we have significantly 
improved our handling on the bulge general properties such as 
distance spread and interstellar extinction, possible variations 
across different regions might still be present. Therefore, it is 
important to quantify the effect of such variations on our mean 
abundance determinations. 

5.1. The bulge distance spread 

Bulge stars are often assumed to be w ell represented by a pop - 
ulation at ~8 kpc from the Sun (e.g. |Rattenburv et al. 2007). 
However, for a wide field study such as the one presented here, 
we must include two factors into our analysis: the bar inclination 
angle and the distance spread of stars. As explained in Section 3, 
we have used the red clump method to derive the mean dis- 
tance of the stars towards each line of sight, therefore the in- 
clination angle of the bulge has been fully taken into considera- 
tion. However, by assuming that all stars are located at the same 
distance, corresponding to the distance modulus derived using 

1 http://mill.astro.puc.cl/BEAM/calculator.php 
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the red clump magnitude, we are introducing an error which is 
directly related to the magnitude width of the red clump in the 
luminosity function. The width of the RC luminosity function, 
for a population located at the same distance, will be determined 
by the intrinsic width given by properties such as metallicity 
and age. However, the dominant factor in the case of the bulge 
RC magnitude distribution width will be the distance spread. We 
have therefore used a boothstrapping method in order to estimate 
the uncertainty arising from adopting a single mean magnitude 
and neglecting the width of the bulge RC. We re-calculated the 
mean metallicity after varying the mean magnitude of the red 
clump by a random value, within 2cr of the mean RC K s magni- 
tude. In each case, cr is obtained from the gaussian fit to the lu- 
minosity function in Equation 1 . We repeated this procedure 600 
times for each VVV tile, since one tile corresponds to our res- 
olution element used to derive the mean red clump magnitude. 
Figure [3] shows the dispersion (<Tdi S ) in the mean metallicity ob- 
tained for each VVV tile after the iteration procedure. 




5 0-5 
Galactic Longitude (1) 



Fig. 3. The metallicity dispersion cr rf/s of each line of sight given 
by the individual VVV tiles, after varying the mean K s magni- 
tude for 600 random values within 2cr of the RC Gaussian fit. 



5.2. Residual differential extinction 

Another source of uncertainties in our derivation of the mean 
metallicity is due to a possible presence of residual differen- 
tial extinction. Although a great improvement is provided by 
the hi gh-resolution extinction map presented in iGonzalez et al.l 
(l2012h . the resolution of this map is limited by requiring a suffi- 
cient number of stars to fit the RC color distribution. The corre- 
sponding resolution of the map, under this limitation, was given 
by 2', 4' and 6' arcmin depending on latitude. Variations on 
scales smaller than these resolution elements are still possible, 
which then produce an increase in the width of the RC color 
and, therefore also in the RGB region that we use to derive the 
metallicities. As the reddening becomes larger towards the inner 
regions of the bulge, the residual reddening follows the same be- 
havior. This creates an excess of redder stars in the CMD which 
would end up in a derivation of higher photometric metallici- 
ties for those stars and creating a metallicity gradient as the one 
we observe in our maps. Thus, we must ensure that the gradient 
we observe is not a product of differential reddening but a real 
feature. 



We therefore evaluate the effect of residual differential ex- 
tinction by studying the width of the {J - K s )o color dis- 
tribution of the RC, measured using a Gaussian fit to the 
RC color in each resolution element of the metallicity map. 
The dispersion, crg-K s ) of the Gaussian fit to RC color 
distribution will depend primarily on metallicity dispersion 
(o~(j-K,),Fe) as well as on the presence of differential redden- 
in g (o-(j-K s \red)- We can thus para metrize o- ( j-k s ) as ct ( j-k s ) = 

V^cr-^o + rfj-We + °\j-K s ),reV where ^Q-k,),0 is the uni- 
form intrisic width of the bulge RC population. In the context 
of this test we can then adopt the hypothesis of having neither 
a metallicity nor a metallicity dispersion gradient in the bulge, 
where both <T(j-k,),o and <T(j-K s ),Fe can be assumed to be con- 
stant across the different regions and the single effect of residual 
differential extinction can be investigated. Then, under this as- 
sumption, the squared root difference between the o-(j_ Ks ) value 
in any region and that of the outermost region of our map, where 
the effects of extinction and certainly residual reddening are not 
present, would be a direct measurment of crg-K s \red- We calcu- 
lated the o-(j-K,),red values for each field using VVV field b217 
as our reddening free control field. 

Finally, we applied a Icr^j-K^.red variation to the color of the 
RGB stars of our control field in the outer bulge and re-derived 
the metallicity distribution for the o~(j-K s \ r ed corresponding to 
each field. Figure [4] shows the map of cr Fe rec i given as the differ- 
ence between the original [Fe/H] of the control field and the new 
value, affected by each <T(j-K t ),red- We observe that the residual 
differential reddening does not mask the shape of the metallic- 
ity gradient and, most importantly, the error values in metallicity 
due to differential reddening are almost neglible for most regions 
(less than 0.05). 




Galactic Longitude (1) 



Fig. 4. Error in the mean metallicity derivation from residual dif- 
ferental reddening. Error values are given by the difference be- 
tween the original [Fe/H] of the control field (b217) and the re- 
calculated metallicity after applying the (T(j-x s ),red correspond- 
ing to each field. 



5.3. Extinction law variations 

The reddening E(J - K s ) from the map presented in Paper II 
can be converted to the corresponding extinctions Aj, Ah, and 
Ak s by adopting a corresponding extinction law. Most of the 
studies towards the Galactic bulge ado pt standard ex t inctio n 
laws consistent with a Ry ~3.1 such as Cardel lTet al.l (Q989). 
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However, there are evidences for a non-standard ex tinction law, 
both in the optical and near-IR, towards the bulge. Gould et al. 
(2001;) and lUdalskH d2003l) demonstrated that the optical redden- 
ing law toward the inner Galaxy is described by smaller total- 
to-selective ratios than th e standard val ues measured for the lo- 
cal interstellar medium. lUdalskil (|2003) measured a value of 
dAj — dE(V -I) = 1.1 towards several bulge fields using OGLE 
photometry, much smaller than the dA; = dE(V - I) = 1.4 
sugg ested by the standard interstellar extin ction curve of Ry= 
3:1 dCardelli et al.1 119891: lOTtonnelll 11994 . Furthermore, not 
only was the optical reddening law toward the bulge found to 
be non-standard, it was also foun d to be rapidly varying be- 
tween sightlines dNataf et al.ll2012l) . The steeper extinction law 
toward the inner Galaxy ha s been subsequently co nfirmed with 
observations in the optical dRevnivtsev et al.ll2.Q10h, by an alysis 
of RR Lyrae stars i n OGLE- HI (Pietrukowicz e t al.ll2.012h. and 
also in the near-IR dMshivama et alj|2009t iGosling et al.ll2009t 
ISchodel et alEoiOh . 

An error in our metallicity calculations can be produced 
when adopting a given extinction law for the complete bulge re- 
gion when such a law might no longer be valid in all regions. In 
particular, if a difference in the extinction law is present between 
the outer and inner bulge, a metallicity gradient could also be a 
consequence of this systematic error instead of a real feature. 
Therefore, we have derived our metallicity map using the stan- 
dard law of Cardelli et al. (1989), shown in Fig. [2] as well as us- 
ing that of Nish ivama et al. (2009) which is now shown in Fig. [5] 
The difference between the metallicity obtained using both red- 
dening laws allow us to evaluate its effect on the observed gradi- 
ent. Figure [6] s hows the diff e rence between th e metallicity map 
derive d using ICardelli et aD (fl989) and that of Nishivama et al.l 
(2009). The maximum difference is observed to be of 0.10 dex 
in the mean metallicity and, as expected, it correlates with the 
amount of extinction. The variations in extinction law are there- 
fore not enough to produce the observed gradient and it also does 
not follow the same spatial behavior of the metallicity gradient. 
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. Metallicity map for the Galactic bulge, as shown in Fig [2] 
the extinction law from Nishiyama et al. (2009) 
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Fig. 6. Difference, [Fe/H] nishiyama - |~Fe/H"| o; m Mii, between the 
metallicities deriv ed using the extinctio n law oflNishivam a et al.1 
d2009h and that of ICardelli etail (1 19891) for the Galactic bulge. 



hand, since the amount of reddening is not necessarily symmet- 
ric with respect to the plane, and it is in fact observed to be higher 
at positive latitudes t han at negative ones dMarshall et al.l i2006: 
iGonzalez et al.ll2012l) . an asymmetric metallicity map would re- 
sult if reddening and metallicity are not being properly disenta- 
gled because of, for example, the use of a non-suitable extinction 
law. Figure Q shows the mean difference found between metal- 
licities above (b > 0) and below (b < 0) the Galactic plane as a 
function of distance from the plane. Error bars show the disper- 
sion in this difference in metallicity for each latitude bin. 




-0.1 
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Fig. 7. Difference in the mean metallicity at different latitudes 
above ([Fe/H] + ) and below ([Fe/H]_) the Galactic plane. The 
me an difference is show n as a doted line. Results obtained using 
the|Cardelli et al. ( 1989) extinction law are shown in black and 
in red for those where the Nishiyama et al. (2009) law is adopted. 



Additionally, we can evaluate the possible effect of extinc- 
tions law variations, in our derived photometric metallicities. 
For this we must first assume that the metallicity distribution 
function is symmetric with respect to the Galactic plane. If this 
is the case, then the difference between metallicities measured 
at positive and negative latitudes should be zero. On the other 



We found a small systematic shift from zero (A[Fe/H] ~ 
-0.01), in the mea n metallicity difference, when using the 
Carde lli et al. I (fl989l) extinction law, which is no longer present 
when the Nishiyama et al. (2009) law is adopted. This particu- 
lar analysis would be independent of the adopted extinction law 
only if the amount of reddening was the same at both positive 
and negative latitudes. Therefore, we can interpret this result as 
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the combined effect of the large amount of reddening at posi- 
tive latitudes and the Cardelli et al. law being less adequate to 
correct for extinction effects in these inner Bulge regions. This 
effect is already included in our previously estimated errors for 
each resolution element of our metallicity map but we can now 
confirm the actual improvement of using the Nishiyama et al. 
(2009) extinction law in the inner Bulge regions (b < 5). 

5.4. The metallicity error map 

After evaluating the individual sources of errors from the dis- 
tance uncertainty, differential reddening, and the extinction law, 
we can obtain a map with the corresponding error for each of the 
subfields in our metallicity map, calculated as: 



cr Fe 



4 



®~Fe,dist ® Fe,red ®~Fe,law 



(4) 



Figure [8] shows the error metallicity map for each subfield. 
This error does not include the method error associated to obtain 
photometric metallicities. However, although the individual er- 
rors of photometric metallicities are expected to be larger than 
those from spectroscopy, we see differences lower than 0.05 dex 
between high resolution spectroscopic and photometric metal- 
licities, when the mean value of the metallicity distribution is 
considered. We therefore consider this method error in the mean 
to be neglictible compared to those previously addressed. 
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Fig. 8. Error map for the mean photometric metallicity for each 
subfield. 



6. The formation scenario of the bulge 

In the recent years, several authors have suggested that the shape 
of the metallicity distribution, and its variation across the bulge, 
is determined by the different contribution of two (or more) com- 
ponents. However, the origin of these components, as well as 
their main characteristics such as their age and metallicity, are 
not yet clear. Therefore, the complete view of the metallicity gra- 
dients accross the bulge is of high importance for constraining 
Galactic models. 

The presence of a radial metallicity gradient is commonly 
considered as a strong argument in favor of the classical bulge 
scenario. In our complete metallicity map we observe such gra- 
dient: a [Fe/H]~0.28 dex/kpc, value that is in remarkable agree- 
ment with the predicted gradient for bulge formed via dissi- 



pational collapse i n the chemo-dynamical model discussed in 
Grieco et al.l(l2012l) . 

However, structural and kinematical studies hint towards a 
different scenario. There is sufficient evidence for the Milky 
Way bulge structure to be consistent with a boxy-bulge, 
formed as a consequence of the vertical buckling of the bar. 
Several photometric studies have shown how th e bar can be 
traced at latitudes well above t he Ga l actic plane (|Stanek et ail 



19971: lLopez-Corredoira et al.l I2005L iRattenburv et al.l 12007 



Dwek et all 1995b iBinnev et al.l 1 1997h . In plane studies of the 



bar, in particular those based on near-IR star counts, have sug- 
gested the presence of a longer bar, in the plane of the Milky 
Way, which extends to larger longitude s and has a different 



orientation angle tha n the o ff -plane bar (Beniamin et al 
Cabrera-Lavers et al.l 120071 12008b IChurchwell etal 



2005; 
20091) . 



Martinez-Valpuesta & Gerhard (2011) have recently presented 
a dynamical model, where a bar and boxy-bulge are formed 
through secular evolution, conciliating both long and thick 
bar components observed in the Milky Way into a single 
original structure. Additiona lly, the observed cylindrical rota- 
tion of this structure led IShen et aTl (120101) to suggest that 
no additional classical component is required to reproduce 
the Milky Way bulge kine mati cal properties. Fu rthermore, 
iMcWilliam & Zoccalil (120 10) and iNataf et all (l2010h presented 
the observations of a double RC feature in the luminosity func- 
tion of the bulge and i nterpreted it as two concentrat ion of stars 
at different distances. iMc William & Zoccalil (1201 Ol) suggested 
that this is the result of the bulge being X-shaped. This finding 
was later followe d by a more com plete mapping of the double 
RC in the work of lSaito et al.l d201 ll) . who provided a more gen- 
eral view of the X-shape morphology of the bulge. 

These conclusions regarding the morphology and kine- 
matics of the Galactic bulge are consistent with results 
from simulations where a bar suffering from buckling in- 
stabilities results in the forma ti on of a peanut or X- s haped 
structure dAthanassoulal 12005b iDebattista & Williamsl [2004; 
iDebattista et alJl2006HMarfinez-Valpuesta et al.ll2006l) that puffs 
up from the plane of the galaxy. When these vertically thick- 
ened bars are observed with an int ermediate angle orientation , 
they would appear as boxy-bulges (jK uijken & Merrifield 1995; 
iBureau & Freemanlll999l:lPatsis et al.ll2002HBureau et al.ll2006l) . 
Additionally, we see from the metallicity map in Fig.|2]that stars 
at larger latitudes are more metal-rich at positive longi tudes with 
respe ct to negative ones. Following the discussion in ( Saito et al. 
1201 ll) . at a given latitude, stars at negative longitudes are further 
away from us, because of the bar viewing angle, and therefore 
their true vertical distance from the plane will be larger than for 
stars at the near side of the bar. The observed horizontal pattern 
in the metallicity map of Fig. 2 would then be fully explained by 
the vertical metallicity gradient along the boxy-bulge. 

In order to reconcile these chemical and structural proper- 
ties of the Galactic bulge, we first need to understand whether 
boxy-bulges alone can al so show a metallicity grad ient as the 
one observed in our map. iBekki & Tsuiimotol (1201 ll) discussed 
this possibility, based on a Galactic model where a boxy-bulge 
is formed out of a two component disk. In their simulation, they 
are able to reproduce a vertical metallicity gradient when the 
bulge is secularly born from an original disk with its own ver- 
tical metallicity gradient, i.e. a metal-poor thick and metal-rich 
thin disk component, such as those observed in the Milky Way. 
The initial chemical properties of the stars in the disk would 
then produce the observed metallicity gradient, as the metal-poor 
stars that dominate the dinamically hotter thick disk would not 
be influenced by the mixing within the bar. 
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Th e boxy-bulge model of Martinez-Valpu esta & Gerhardl 
(1201 ll) . which has shown to succesfully reproduce the stellar 
count profiles for both the inner and outer bulge, also shows 
a remarkable agreement with the complete metallicity gradient 
shown in this work (Martinez-Valpuesta et al. 2013, in prepa- 
ration). We recall that the model does not include a classical 
component, hence demonstrating that the buckling process of a 
secularly evolved bar could be a sufficient process to explain the 
formation history of the bulge. 

In support of the above scenario, there is the chemical 
similarity observed between bulge and disk stars. Metal-rich 
stars of the bulge are observed to be as alpha-poor as those 
of the thin disk, while metal-poor bulge stars show an alpha- 
element enhancement which is indistinguishable from that of 
the Galactic thick disk and dominate at high e r scale heights 
dMelendez et all l2008t lAlves-Brito et alJ l2010t iGonzalez etall 
201 la] ). Furthermore, the mainly old stellar age s observed in the 
bulge (IZoccali et al.l2003"HcTa rkson et al. 201 1) place additional 
constrains to such formation scenario. The first original thick 
disk would have to be necessarily formed fast (~1 Gyr), most 
likely as the result of dissipative collapse or the merging of gas- 
rich subgalactic clumps, in order to be consistent with the early 
formation of the hotter, old and alpha-enhanced component. The 
inner stellar disk would then need to proceed with its star for- 
mation, relatively fast compared to the outer thin disk, to then 
become the inner parts of the boxy-bulge. 

Finally, we recall that recent studies have also suggested that 
the chemodynamical properties of the bulge are consistent with 
a change of the dominating population between the outer and 
inner bulge: a metal-poor, alpha-rich component with spheroid 
kinematics which is homogeneous along the minor axis, and 
a meta l-rich, alpha-poor bar population domi n ating the inner 
bulge dBabusiaux et all l20~iot iHill et alJ l20Tlt IGonzalez et alJ 
1201 lah . Under this interpretation the vertical metallicity gradi- 
ent will be naturally produced by the density distribution of both 
components at different lines of sight. Although this is also con- 
sistent with the general view seen in our complete metallicity 
map, the large individual metallicity errors in our photometric 
method in comparison to those from spectroscopy, do not allow 
to invenstigate further the bimodality in the metallicity distribu- 
tion (see Fig. 1). Therefore, the presence of a bimodal metallic- 
ity distribution, and its correlation with kinematics, remains to 
be spectroscopically confirmed within a large spatial coverage. 

7. Conclusions 

We used the RC stars of the Galactic bulge as a standard candle 
to build a reddening map and determine distances. This allowed 
us to construct the color magnitude diagrams of bulge stars in 
the absolute plane (Mg s ,(J - K s )o) and then to obtain individ- 
ual photometric metallicities for RGB stars in the bulge through 
interpolation on th e RGB ridge lines of cluster sample from 
IValenti et alJ J2007), taking into account the population correc- 
tions for the RC absolute magnitude as a function of metallicity, 
in an iterative procedure. The final metallicity distributions were 
compared to those obtained in the spectroscopic studies show- 
ing a remarkable agreement. The main result is the first high 
resolution metallicity map of the Milky Way bulge with a global 
coverage from -10 < / < +10 and -10 < b < +5, excluding 
the central area {\b\ < 2) where too high crowding and strong 
differential extinction compromise the measurements. 

The global bulge metallicity map is consistent with a verti- 
cal metallicity gradient (~ 0.04 dex/deg), with metal-rich stars 
([Fe/H]~0) dominating the inner bulge in regions closer to the 



galactic plane (\b\ < 5). At larger scale heights, the mean metal- 
licity of the bulge population becomes more metal-poor. This fits 
in the scenario of older, metal-poor stars dominating the bulge at 
higher scale heights above the plane as a consequence of original 
gradients present in the Galactic disk. Therefore our metallicity 
map supports the formation of the b oxy bulge from the buck- 
ling instabillity of the Galactic bar (iBekki & Tsuiimotoll201 lb 
Martinez-Valp uesta & Gerhardll201 ll ). This scenario is also sup- 
ported by the chemical similarities observed between the metal- 
poor bulge and the thick disk stars. 
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